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Quantum systems can provide outstanding performance in various sensing applications, ranging from 
bioscience to nanotechnology. Atomic-scale defects in silicon carbide are very attractive in this respect 
because of the technological advantages of this material and favorable optical and radio frequency 
spectral ranges to control these defects. We identified several, separately addressable spin-3/2 centers in 
the same silicon carbide crystal, which are immune to nonaxial strain fluctuations. Some of them are 
characterized by nearly temperature independent axial crystal fields, making these centers very 
attractive for vector magnetometry. Contrarily, the zero -field splitting of another center exhibits a giant 
thermal shift of —1.1 MHz/K at room temperature, which can be used for thermometry applications. 
We also discuss a synchronized composite clock exploiting spin centers with different thermal 
response. 



Nanoscale magnetic field and temperature sensing based on quantum properties of the nitrogen-vacancy 
(NV) defect in diamond 1 has been proposed 2 " 5 and demonstrated at ambient conditions 6 " 8 . This sensing 
technique can potentially allow for monitoring neuron activity 9 , imaging of single proton spins in 
complex molecular systems 1011 and measuring the heat produced by chemical reactions inside living cells 12,13 . 
The principle of sensing is based on the optical readout of the spin resonance frequency subject to the crystal field 
D and magnetic field B. 

These successful experiments have stimulated the search for other solid-state spin systems with similar abilities. 
Atomic-scale defects in silicon carbide (SiC) are attractive in this respect. Apart from the obvious technological 
opportunities due to the well-developed device fabrication protocols in SiC, these defects reveal favourable optical 
and radio frequency (RF) spectral ranges for their control. Namely, (i) they are optically active in the near infrared, 
characterised by a deep tissue penetration and a weak absorption in optical glass fibers. And (ii) they can be 
controlled by the MHz RF-fields, used in standard magnetic resonance imaging systems. Because of about 250 SiC 
polytypes are known, there should exist more than thousand different spin defects in SiC with distinct character- 
istics 1415 . One can select defects with the most suitable properties for a concrete task, which is not possible for one 
universal sensor. 

Here, we identify separately addressable spin centers within the primary intrinsic defects in 6H-SiC and 4H- 
SiC. Some of these half-integer spin (S = 3/2) centers are insensitive to nonaxial strain as a consequence of the 
Kramers degeneracy theorem and reveal within our accuracy temperature independent crystal field. Thermal/ 
strain fluctuations may limit the magnetometry performance. In particular, thermal fluctuations lead to fluctua- 
tions of the spin resonance frequency because D is generally a function of temperature (T). The NV defect 
demonstrates the thermal shift of the spin resonance frequency ft = — 74 kHz/K 16 , and temperature fluctuations 
of 1°C results in a magnetic field uncertainty of a few /nT. By stabilizing temperature and using advanced 
measurement protocols in a bias magnetic field, this uncertainty can be significantly reduced 17 , while in case 
of SiC this is already the intrinsic property. 

On the other hand, a large thermal shift is necessary for temperature sensing. We also identify spin-3/2 centers, 
demonstrating giant thermal shifts up to = — 1.1 MHz/K at ambient conditions. This is 14 times larger than 
that for the NV defect and hence can potentially be used to enhance the thermometry sensitivity. 
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Figure 1 | Optical detection of spin defects in 6H-SiC. (a) Crystallographic structures of the Schottky [V si -V c , V si (V2), V si (V3)] and Frenkel (V si -Sii) 
defects. All the defects have the symmetry axis oriented along the c-axis of the crystal lattice, (b) Room-temperature zero-field ODMR spectrum. 
Inset: optical pumping scheme of spin defects. 



I. Room-temperature ODMR of intrinsic defects in 
6H-SiC 

The Schottky and Frenkel defects are the primary stoichiometric 
defects in solids. A divacancy (V Si -V c ) -consisting of chemically 
bound silicon vacancy (V S i) and carbon vacancy (V c )-is an example 
of a Schottky defect in SiC [Fig. 1(a)] being extensively investi- 
gated 18 " 21 . Another example of a Schottky defect is an isolated silicon 
vacancy 18,22 " 28 perturbed by the nearest carbon vacancy along the c- 
axis of the SiC lattice 15 . There are several possible configurations and 
two of them are labeled in Fig. 1(a) as V Si (V2) and V Si (V3). Frenkel 
defects have also been identified in SiC and one of such defects-a 
silicon vacancy and a corresponding interstitial Si atom located at a 
distance of 6.5 A along the c-axis (V Si -Sii) 29 -is shown in Fig. 1(a). 

A common feature of these defects in SiC is that they have a high- 
spin ground state 19,29 . The corresponding spin Hamiltonian is written 
in the form 

n = g e fi B BS + D(T)[S 2 z -S(S+l)/3]. (1) 

Here, g e ~ 2.0 is the electron g-factor, \i B is the Bohr magneton and S z 
is the projection of the total spin S on the symmetry axis of the defect 
(in case of the defects considered here it coincides with the c-axis of 
SiC). Without external magnetic field (B = 0) the ground state is split 
due to the crystal field D and for S = 3/2 the zero-field splitting (ZFS) 
between the m s = ±1/2 and m s = ±3/2 sublevels is equal to 2D [the 
inset of Fig. 1(b)]. Remarkably, the zero-field splitting is an indi- 
vidual fingerprint for each defect. In table I we summarize the ZFS 
values for different intrinsic defects in the most commonly encoun- 
tered polytype 6H-SiC. 

We now demonstrate, using optically detected magnetic res- 
onance (ODMR) technique, that defect spins in 6H-SiC can be initi- 
alized and read out at ambient conditions, which is the basis for 



various sensing applications. Optical excitation with 785 nm fol- 
lowed by spin- dependent recombination through the metastable 
state (MS) results in a spin polarization of the defect ground state 
[the inset of Fig. 1(b)]. On the other hand, photoluminescence (PL) 
rates from the excited state (ES) to the ground state (GS) involving 
different spin sublevels are different. When a resonance RF equal to 
the ZFS is applied (v 0 = 2D/h for S = 3/2), it induces magnetic dipole 
transitions between the spin-split sublevels resulting in a change of 
the photoluminescence intensity (APL). Since single photons emit- 
ted by a single defect can be detected, this ODMR technique is now a 
standard method to probe single NV spins in diamond 1 , which has 
the resonance frequency around v 0 = 2.87 GHz. 

In order to probe spin defects in SiC we have extended the RF 
spectral range from a few GHz down to a few tens of MHz, usually 
assigned for nuclear magnetic resonance experiments. Figure 1(b) 
shows a typical ODMR spectrum, i.e., relative change of the photo- 
luminescence intensity APL/PL as a function of applied RF. Two spin 
resonances at v 0 = 28 MHz and v 0 = 128 MHz agree well with ZFS 
of the silicon vacancy defects V Si (V3) and V Si (V2), respectively 
(table I). Another spin resonance is observed at v 0 = 367 MHz 
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Figure 2 | Vector magnetometry. (a) Zeeman splitting of the ±3/2 (upper) and ± 1/2 (lower) spin sublevels as a function of the angle 6 between magnetic 
field B = 0.1 mT and the defect symmetry axis (2D = 128 MHz). The arrows indicate four possible RF-induced transitions, (b) Comparison of the 
relative shifts of the spin resonance frequencies dv/dB as a function of 6 for S = 1 and S = 3/2. 



and we ascribe it to the Frenkel pair V Si -Sii (table I). A discrepancy of 
10% can be explained by the strong dependence of its ZFS on tem- 
perature and laser power, as described later in this work (see also 
Supplemental Material). Finally, we discuss the origin of the spin 
resonance at v 0 = 1.03 GHz. We find that this spin resonance is 
strongly temperature dependent-D(T) is discussed later in detail- 
and at T = 10 K we measure the zero -field splitting to be v 0 = 
1.22 GHz. This agrees reasonably well with ZFS of the divacancy 
VsrV c (table I). On the other hand, we register PL up to 1050 nm, 
but the V Si -V c PL has maximum above 1 100 nm 20 . We would like to 
mention that the magnetic field dependencies (presented later) un- 
ambiguously indicate S = 3/2 for this defect, while in the earlier 
studies the S = 1 ground state of the divacancy was assumed 19 . For 
this reason we label it as unidentified defect (UD) in all figures, 
indicating that its origin should be carefully examined. 

First, as a characteristic example, we consider theoretically the 
Zeeman splitting of the V Si (V2) spin sublevels for different orienta- 
tions (given by the polar angle 0) of a weak magnetic field 
B<2D/g e ix B [Fig. 2(a)]. The V Si (V2) defect has spin S = 3/2 ground 
state 15 ' 22 ' 24 and the upper m s = ±3/2 spin sublevels are split as 

±-g e ft B B cos0. Interestingly, the lower m s = ±1/2 spin sublevels 
are mixed due to the perpendicular-to-the-symmetry-axis field com- 
ponent and split differently as ±^g e fi B B(l + 3 sin 2 O) 1 ^ 1 , which 

results in four possible RF-induced transitions. We obtain the prob- 
ability of these transitions [color-coded in Fig. 2(b)] as 

^H(j|BiS|fc>| 2 , (2) 

where \j), \k) are the eigenstates of (1) and B Y ±c is the driving RF 
field. 

Figure 2(b) compares the relative shifts of the spin resonance 
frequency dv/dB depending on the magnetic field orientation for 
S = I and S = 3/2. In case of a spin-1 system (i.e., the NV defect) 
it follows \g e fiB cos 6/h\ (h is the Planck constant) and tends to zero 
for B±c (0 = 90°). Therefore, using a single (or equally oriented) 
spin-1 defect(s) the field projection on the defect symmetry axis B z = 
B cos 6 rather than the absolute value is detected. In fact, due to the 
presence of the nuclear spin bath, the NV spin-echo decay curve 
contains characteristic "collapses and revivals" from which the total 
magnetic field can be determined 30 . However, the nuclear gyromag- 
netic ratio is much smaller (by a factor of 0.0003 for 13 C) than that of 
electron, making this method less sensitive. Furthermore, this 
method is not applicable for isotopically purified crystals, providing 
the highest sensitivity. In contrast, in case of a spin-3/2 system there 
are two pairs of resonances. Therefore, from the comparison of the 



relative ODMR frequencies and amplitudes both the absolute value B 
and its orientation with respect to the defect symmetry axis can be 
reconstructed. 

In order to demonstrate this advantage of the spin-3/2 defects 
in SiC for magnetic sensing applications it is necessary to ascertain 
their symmetry axes. This can be done in high magnetic fields 
B^>2D/g e /n B . In this case, a pair of spin resonances should be 
observed for each spin defect, the separation between those being a 
function of 9. We have performed this experiment in the X-Band at a 
fixed frequency v = 9.4 GHz while sweeping the magnetic field and 
an ODMR spectrum obtained for B\ \c is shown in Fig. 3(a). Indeed, 
four pairs of spin resonances from four distinct spin defects are 
observed, in accord with four lines in the zero -field ODMR spectrum 
of Fig. 1(b). The resonance magnetic fields B ± related to each pair are 
presented in Fig. 3(b) as a function of the angle between the field 
direction and the c-axis. 

According to spin Hamiltonian (1), the difference between these 
fields is given by 

AB = B + -B_ = 2 RiIl (3cos 2 0-l). (3) 

As Eq (3) describes the experimental data quite well [Fig. 3(b)], one 
immediately concludes that all the defects under consideration are in 
line with the c-axis, as also depicted in Fig. 1(a). Furthermore, for B\ \c 
(0 = 0°) ZFS can independently be measured 27 as 2D = g e fi B AB/2 
and from the comparison with the zero-field data of Fig. 1(b) all the 
resonances in Fig. 3 are unambiguously assigned. 

Having established the orientation of the spin defects in our 6H- 
SiC sample, we apply a weak magnetic field B = 2.2 mT perpendic- 
ular to the defect symmetry axis (6 = 90°) [Fig. 4(a)]. In this con- 
figuration the difference between different spin systems is most 
pronounced [see Fig. 2(b)]. For a spin-1 defect such oriented weak 
magnetic fields should have no effect. In contrast, we observed large 
Zeeman splittings for all four ODMR lines [symbols Fig. 4(b)]. The 
experimental data are perfectly described by the spin Hamiltonian 
(1) with S = 3/2 and corresponding transition rates of Eq. (2), as 
shown by the solid lines in Fig. 4(b). For comparison, the calculations 
with S = 1 are shown by the dotted lines. Remarkably, for V Si (V2) 
and especially for V Si (V3) the magnetic field behavior is not linear 
because the Zeeman splitting in a magnetic field of a few mT is 
comparable or larger than ZFS in these silicon vacancy defects, which 
results in the intermixing of all four spin sublevels. According to the 
magnetic field dependencies of Fig. 4(b), all defects under considera- 
tion have the S = 3/2 ground state. This is in agreement with the 
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Figure 3 | Axial symmetry of spin defects in 6H-SiC. (a) Room-temperature ODMR spectrum in the X-band (v = 9.4 GHz) as a function of 
magnetic field BWc. (b) ODMR frequencies for different angles between the magnetic field B and the oaxis of 6H-SiC. Experimental data are shown by 
symbols, the solid lines correspond to the calculations to Eq (3). 



earlier ascertained spin structure of the V Si defects 15 and with the 
model of the V Si -Sii Frenkel pairs 29 . 

We now discuss the effect of strain. In case of S = 1 spin system, as 
the NV defect in diamond and the divacancy in SiC, the local strain 
lifts up the degeneracy of the m s = ± 1 spin sublevels even without 



external magnetic field. It is described by an additional term 
E^S 2 x — S^j in the spin Hamiltonian, where E is the transverse ZFS 

parameter. This may limit the magnetometry sensitivity because the 
spin -splitting in external magnetic field for E 0 is described 



V/- 




Frequency, v (MHz) 

Figure 4 | Effect of external magnetic field in 6H-SiC. (a) Room-temperature ODMR spectrum obtained in a magnetic field B = 2.2 mT. (b) ODMR 
frequencies as a function of magnetic field. Experimental data are shown by symbols, the solid and dotted lines correspond to the calculations for 
S = 3/2 and S = 1, respectively. In all panels B _L c. 
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Figure 5 | Temperature effect in 6H-SiC. (a) ODMR spectra in the X-band (v = 9.4 GHz) recorded at different temperatures, (b) Zero-field 
splitting 2D for different spin defects as a function of temperature. The solid lines are third order polynomial fits. Inset: The same as (a), but rescaled 
to highlight a strong UD temperature dependence. 



by y (g e /j, B B z ) + £ 2 , i.e., the relative shift of the spin resonance 
frequency is quadratic rather than linear for small B. By improving 
the quality of bulk diamond material, the parameter E can be reduced 
below 1 MHz and further eliminated in a small bias magnetic field. 
However, magnetic fields may lead to additional inhomogeneity, 
which is in some experiments unwanted. Furthermore, in nanodia- 
monds, which are very interesting for bio-sensing, strain is usually 
present and the transverse ZFS parameter can be as large as E = 
20 MHz 31 . The situation is qualitatively different for half-integer spin 
systems, like spin defects with S = 3/2. According to the Kramers 
theorem, the double degeneracy of the m s = ±1/2 and m s = ±3/2 
states can only be lifted up by an external magnetic field, meaning 
that essentially £ = 0. This makes these defects (in particular, silicon 
vacancies in SiC) more robust against strain fluctuations. Excitons 
and trions 32 are a well-known example of this consequence of the 
Kramers degeneracy theorem in solid-state system. Excitons (with 
integer angular momentum 7=1) demonstrate a fine structure 
splitting in quantum dots. It is caused by local strain due to lattice 
mismatch between the material of the quantum dot and the barrier. 
On the other hand, the trions (with half- integer angular momentum 
/ = 3/2) reveal no fine structure splitting under the same strain 
conditions 33 . 

In order to examine the effect of temperature fluctuations, we have 
measured ODMR spectra in the temperature range from 10 to 320 K 
[Fig. 5(a)]. The experiment is performed in the X-band and D for 
each defect is determined using Eq. (3). As one can see from Fig. 5(b) 
the parameter D of the V Si (V2) and V Si (V3) defects is temperature 
independent within the accuracy of our experiment (a few kHz/K). 
This is an additional strong argument to use these defects for 
magnetometry. 

In contrast to V S i(V2) and V S i(V3), for the Frenkel pair V S i-Sii we 
observe a reduction of its ZFS of ca. 50%, from 0.6 down to 0.3 GHz 



when the temperature increases from 10 to 320 K [Fig. 5(b)]. The 
corresponding temperature dependence is well described by a third 
order polynomial 2D(T) = a 0 + a{T + a 2 T 2 + a 3 T 3 with a 0 = (667 
± 4) MHz, a x = (-1.7 ± 0.1) MHz/K, a 2 = (2.3 ± 0.8) X 10" 3 
MHz/K 2 , and a 3 = (-8 ± 15) X 10" 7 MHz/K 3 . The defect with the 
largest ZFS (labeled as UD in Fig. 5) has a slightly weaker temper- 
ature dependence (17% relative reduction of D) but reveals the stron- 
gest thermal shift IdDldT at room temperature [the inset of 
Fig. 5(b)]. From the polynomial fit with a 0 = (1222 ± 2) MHz, a x 
= (2.4 ± 0.4) X 10" 1 MHz/K, a 2 = (-5.2 ± 0.3) X 10" 3 MHz/K 2 , 
and a 3 = (6.5 ± 0.5) X 10" 6 MHz/K 3 we find ft = IdDldT = 
-1.1 MHz/K at T = 300 K. This value is 14 times larger than the 
thermal shift of the NV defect in diamond 16 . 

II. Intrinsic defects in 4H-SiC and magnetic field 
sensing 

We also observe very similar behaviour in 4H polytype of SiC, which 
can be grown of high quality. Figure 6(a) shows a room-temperature 
ODMR spectrum of 4H-SiC, obtained in the X-Band spectrometer. 
Following the angle dependencies presented in Fig. 6(b) we ascribe 
the ODMR lines to the silicon vacancy V Si (V2) (with ZFS of 
67 MHz) and the Frenkel pair V Si -Sii (with ZFS of 408 MHz) 23 ' 27 ' 29 . 

We now concentrate on the ODMR spectra of the V S i(V2) defect in 
zero and weak bias magnetic fields [Fig. 7(a)]. In zero magnetic field, 
the ODMR line has maximum around v 0 ~ 70 MHz [Fig. 7(b)]. Upon 
application of a magnetic field along the c-axis (6 = 0°), the ODMR 
line is split and one observes two resonances [Fig. 7(a)]. They shift 
linear with magnetic field as v B = |v 0 ± g e ^BB/h\ [Fig. 7(b)], in accord 
to Hamiltonian (1) for B||c. From the linear fit we find the V S i(V2) 
ZFS in our crystal with higher accuracy v 0 = 70.2 ±0.3 MHz. 

In our earlier studies we have reported that the ODMR linewidth 
in SiC strongly depends on the RF power 15 . Therefore, to achieve the 
narrowest linewidth we use the lowest RF power, which still provides 



SCIENTIFIC REPORTS | 4:5303 | DOI: 1 0.1 038/srep05303 



5 



CL 

< 

CL 







4H-SiC 






T=305K 


V -Si. 

Si i 






♦ J 




V * 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


V 



-40 
-20 

C ok 

^ 20 
< 40- 
60- 
80 



2x408 MHz 



2x67 MHz 




290 300 310 320 330 340 350 
Magnetic field, 6 (mT) 



360 370 380 



Figure 6 | Room-temperature ODMR in 4H-SiC. (a) An ODMR spectrum in the X-band (v = 9.4 GHz) as a function of magnetic field Bile 
(b) ODMR frequencies for different angles between the magnetic field B and the oaxis. Experimental data are shown by symbols, the solid lines 
correspond to the calculations to Eq (3). 



a reasonable signal-to-noise ratio. Additionally, we apply a small 
bias magnetic field of 3.5 mT to eliminate the broadening due to 
environment magnetic fields. The result is presented in Fig. 7(c) 
and can be well fitted to 



APL 



r„ 



[2nT*{v B + nA Si /2-v)Y + l 



(4) 



Here, A Si is the hyperfine constant describing the interaction 
ft = IAS with the next-nearest-neighbor (NNN) 29 Si (I = 1/2) nuc- 
lear spins. We find A Si = 9.1 ± 0.2 MHz, which is comparble with 
the earlier reported value 23 . The relative amplitudes of the peaks 
r 0 : Tj : r 2 is given by the probability to find n spin-active 29 Si iso- 
topes among the 12 Si NNN. We find 0.89 X 10" 4 : 0.24 X 10" 4 : 0.36 
X 10~ 5 , which are close to the ratios expected for the natural abund- 
ance of 29 Si (4.7%). 

The resonance linewidth in Fig. 7(c) is inversely proportional to 
the inhomogeneous spin coherence time, which is T\ = 120 ns 
according fit to Eq. (4). The accuracy to measure the relative change 
of the frequency (Sv) is given by the maximum value of the derivative 
(<5r '/ 'Sv) max = tiT^Tq. With the root-mean-square deviation c>r = 
3.5 X 10" 6 [as in Fig. 7(c)] we obtain 3v = 100 kHz. Taking into 
account the integration time of 8 s per point, this corresponds to the 

magnetic field sensitivity 3B = 10fiT ^ a/Hz. The sensitivity can be 

improved by increasing of T 2 . It is limited by the spin-lattice relaxa- 
tion time, which is T x = 100 fis 26 for V S i in SiC and can be achieved 
using spin-echo techniques. Further improvement by approximately 
an order of magnitude should also be possible in isotopically purified 
28 Si 12 C, as was proposed recently 27 . For instance, the record spin 
coherence time of the NV defect in 12 C-diamond is T 2 = 1.8 ms 34 . 
In order to compare the sensitivity of the V Si defect with the sens- 
itivity of the NV defect, we first estimate the product of the overall 
detection efficiency (C) and the number of the V Si (V2) defects (IV) in 
our experiments, by taking into account the signal level (ca. 10 nW) 
and the PL lifetime (6 ns, ref. 35). We then estimate the expected 
sensitivity for optimized samples to be on the order of magnitude 

(55=10 nT/VcN-Hz. 



Our experiments have been performed on an ensemble of defects. 
Very recently, the detection of some single defect centers in SiC has 
been reported 36 , and we believe that the same approach can be used to 
isolate single V Si -related defects. An important issue in this case is the 
ODMR contrast. It increases with RF power and we observe T 0 = 
0.11% for RF power of 36 dBm. This may lead to some limitations 
compared to single NV defects with larger ODMR contrast (r 0 ~ 
10%), such as a requirement of longer integration time to achieve the 
equal sensitivity with the same other parameters. On the other hand, 
we collect PL from various types of defects and only one of them gives 
rise to the ODMR contrast. This means that T 0 of a single defect 
should be larger than in ensemble experiments. Alternatively, the 
ODMR contrast can be sizeably increased by a proper choice of 
the spectral detection window 14 . Furthermore, in many sensing 
applications it is not necessary to use just one defect. All the silicon 
vacancy related defects considered here are aligned along the same 
axis, and, as intrinsic defects, they can be created of high concentra- 
tion without need for additional doping. 

III. Temperature sensing 

We now discuss the temperature sensitivity, which is proportional to 
(IdD/dT)' 1 3v and therefore it is natural to select a defect with a large 
thermal shift such as the Frenkel pair V Si -Sii [Fig. 5(b)]. A large 
change of ZFS with temperature cannot be explained simply by 
thermal lattice expansion. Assuming that D is due to dipolar coup- 
ling oc r" 3 between the silicon vacancy and the interstitial Si atom (r 
is the distance between them), one should have dDf(DdT) ~ —3a. 
Here, a = 4.0 X 10~ 6 K" 1 is the thermal expansion coefficient of SiC, 
but we observe two orders of magnitude larger value of 2 X 10" 3 K" 1 . 
A possible explanation is that the Si interstitial (SiJ is not rigid in the 
SiC lattice and a small perturbation may result in a significant shift in 
its position relative to V Si , leading to a large change of ZFS. 

The UD defect demonstrates even larger thermal shift. Using the 
same procedure as for the V Si (V2) defect we estimate the temper- 
ature sensitivity 3T ~ 1 Y^j VHz. Remarkably, the ODMR linewidth 

is significantly larger than for V Si (V2) and can be explained by a high 
sensitivity to the local environment. In high-quality samples or for 
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V si ( V2) ODMR frequencies as a function of magnetic field for B 1 1 c. Experimental data are shown by symbols, the solid lines correspond to the calculation 
as explained in the text, (c) An ODMR line of V si (V2) fitted to Eq (4). Arrows indicate the positions of the 29 Si isotope peaks. 



single defects the inhomogeneous broadening should be suppressed, 
leading to homogenous linewidth l/(nT 2 ). In isotopically engineered 
samples and with the use of advanced measurement protocols 81213,17 

the projected temperature sensitivity is 3T= 1 mK j VCN Hz. 

Another potential application could be a synchronized composite 
clock, as was originally proposed for the NV defects in diamond 37 . 
We find the Allan variance (ST/T 0 )/ (inQ^fx) corresponding to the 
resonance in Fig. 7(c) to be 2 x 10 ~ 4 / y/z. By applying a spin-echo 
technique to isotopically purified samples it should be possible to 
significantly increase the quality factor Q = nT 2 v 0 and with T 2 = 
1 ms the Allan variance approaches 10 -8 /Vt. Further improve- 
ment by several orders of magnitude can potentially be achieved 
by increasing the number of defects (signal-to-noise ratio 
Tq/STcc Vn), as this does not require intentional doping with para- 
magnetic impurities (i.e., nitrogen for the NV defect in diamond). 
In this case, the main limiting factor is temperature uncertainty. 
For instance, the temperature stabilization within AT = 0.01 K 
leads to a fractional frequency stability ATp/v 0 < 10~ 6 for the V Si 
defects in SiC. To compensate temperature fluctuations, it has 
been proposed to lock the frequency difference between two 
defects with different thermal shifts and /? 2 ), resulting in a 

frequency uncertainty scaled as J 1 + 2(/? 1 /A/? 12 ) 2 relative to 
the temperature-insensitive case 37 . This idea can be implemented 
within one SiC crystal by using different types of defects. For 
instance, the use of the V Si (V2) defect (Pi) and the UD defect 
(P 2 ) in 6H-S1C yields ft/A^ < 10" 2 . 

Summarizing, intrinsic defects in SiC demonstrate complement- 
ary properties to the NV defect in diamond. Variety of the defect 
characteristics makes them very perspective for various sensing 
applications, such as magnetometry, thermometry and chip-scale 
timekeeping. An important aspect is the integrability of these defects 
within SiC-based electronic devices 28 . Another intriguing possibility 
is to use biocompatibility of SiC nanocrystals to perform nanonscale 
magnetic field and temperature imaging in biological systems, which 
can be done at a fixed RF frequency as well 38 . 



Methods 

We present experimental results for 6H and 4H polytype SiC crystals grown by the 
standard sublimation technique. In order to generate defects our crystals have been 
irradiated with neutrons (5 MeV) to a dose of 10 16 cm" 2 . 

The full-field ODMR experiments are performed in a X-band cavity (Q = 3000, 
v = 9.4 GHz) with direct optical access. The samples are mounted in a liquid helium 
flow cryostat with a possibility to vary temperature from 10 K to slightly above room 
temperature. The temperature sensor is situated in the cryostat about 1 cm below the 
sample. A diode laser operating at 785 nm (100 mW) is used to optically pump all 
types of defects in SiC. The PL excited by this laser is passed through a 950-nm 
longpass filter and detected by a Si photodiode (up to 1050 nm). 

In the zero-field ODMR experiments, the RF radiation provided by a signal gen- 
erator is amplified up to 36 dBm (4 W) and guided to a thin copper wire terminated 
with 50-Q impedance. The laser beam (15 mW) is focused close to the wire using an 
optical objective. The filtered PL signal (10 nW) is coupled to an optical fiber and 
detected by a Si photodiode. Weak magnetic fields are applied using a permanent 
magnet, the field orientation and strength being ascertained using ODMR on the NV 
defects in diamond as a reference. In both experiments we chop the RF radiation 
(typically below 1 kHz) and the output PL signal is locked-in. 
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